Granular dynamics in compaction and stress relaxation 
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Elastic and dissipative properties of granular assemblies under uniaxial compression are stud- 
ied both experimentally and by numerical simulations. Following a novel compaction procedure 
at varying oscillatory pressures, the stress response to a step-strain reveals an exponential relax- 
ation followed by a slow logarithmic decay. Simulations indicate that the latter arises from the 
coupling between damping and collective grain motion predominantly through sliding. We char- 
acterize an analogous "glass transition" for packed grains, below which the system shows aging in 
time-dependent sliding correlation functions. 
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Mechanically agitated granular materials are charac- 
terized by slow relaxation dynamics, arising from the re- 
arrangement of the constituent grains within the volume 
in which they are confined pj . This leads to a slow com- 
paction of the system volume, which follows a logarith- 
mic decay in time, as seen in several experiments and 
predicted by theoretical models of granular compaction 
[3, This property has prompted analogies between 
the physics of athermal granular materials and thermal 
glasses, the theory of which is better understood at the 
fundamental level . The field of granular matter there- 
fore benefits from such parallels, as new ways of inves- 
tigating the system's complex properties are discovered. 
An advantage of using granular materials over glasses is 
that it facilitates a much easier exploration of the rni- 
crostructure through grain-grain interactions. 

Logarithmic relaxation and rate-dependent strength- 
ening have also been observed in compressed granular 
matter although the underlying mechanism is still 
under much debate. The collective rearrangement of the 
grains in the bulk could be responsible for the slow re- 
laxation, as suggested by experiments on slowly sheared 
granular materials by Hartley and Behringer [4j . On the 
other hand, it is also known that aging occurs at the con- 
tacts between the particles which manifests itself as 
a logarithmic increase of the friction coefficient between 
grains as a function of time. It could also be respon- 
sible for the observed slow dynamics, as has been sug- 
ested in experiments by Ovarlez et al. and Nasuno et al. 
Q , which show rate dependence and slow strengthen- 
ing characteristics of aging at the interparticlc contacts, 
respectively. 

The goal of this Letter is to demonstrate the existence 
of slow relaxation in the response of dense granular mat- 
ter to infinitesimal strain perturbations and to elaborate 
on the origin of the dynamics. The experiments reveal a 
very slow stress relaxation under a constant applied dif- 
ferential strain. This behavior is well characterized by 
a two-step relaxation dynamics, analogous to the slow 
relaxation in "glassy" systems p| . 

We investigate this dynamics via computer simula- 
tions, which employ various dissipative processes into 
the system in order to compare their relative effects. 



The results show that the main process responsible for 
the logarithmic stress relaxation is the collective parti- 
cle motion and rearrangements of grains, predominantly 
through sliding. We compute the fraction of sliding par- 
ticles and find that when the damping in the system ex- 
ceeds a critical value, a slow increase in the number of 
sliding particles as a function of time is observed. This 
slow strengthening leads to the logarithmic stress relax- 
ation. Moreover, the system shows the hallmark of glassy 
behavior: aging in the sliding correlation function. 

Experimental arrangement.- We use an INSTRON 
press to measure the mechanical properties of granular 
assemblies confined in a cylindrical cup and piston. The 
machine was strain controlled and enabled oscillatory, 
step, and ramp comprcssional tests up to a maximum 
limit of 300kN load. The cylinder cup had the following 
dimensions: diameter d =5.08 cm, height h =7.62 cm 
and wall thickness I = 0.85 cm. These dimensions were 
chosen to achieve a good statistical ensemble for nearly 
monodisperse glass beads of diameter (355±5)^to. We 
perform two tests: first, we measure the strain response 
(compaction) under large stress oscillations, and then we 
measure the stress response under an infinitesimal step- 
strain perturbation. 

Compaction experiments (stress-controlled).- Before 
performing the stress relaxation experiments it was nec- 
essary to develop a method of reaching the jammed 
state, ensuring that reproducible experiments can be per- 
formed. 

Here we introduce an alternative method to the one of 
Knight et al. of compaction based on oscillatory pres- 
sure of varying amplitude to generate reversible jammed 
states. The compaction procedure consists of the stages 
depicted in Fig. ^ where we show the displacement of 
the piston Ah during the application of the stress a. The 
material is first compressed with a constant, slow veloc- 
ity of 0.1 mm/min to a target stress of 0.5 MPa (from 
A to B in Fig. ^ which shows only the strain response), 
then oscillated between zero and double the mean value 
with a frequency of 1Hz (from B to C, the amplitudes of 
oscillation are always taken as equal to the mean value). 
Each compaction procedure consists of 3,600 cycles (un- 
til C), after which the material is again slowly released 
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FIG. 1: Compaction procedure: Sinusoidal displacement Ah 
of the system height in response to an oscillatory stress around 
the mean value of 0.5 MPa (from B to C) and 1 MPa (from E 
to F). The solid black line is the window average displacement 
of the height and shows the logarithmic behavior of the strain 
under large stress amplitude according to Eq. Q for the later 
stages of the time evolution. The shaded gray area symbolizes 
the amplitude of the oscillatory strain response. The inset 
shows a schematic of the experimental arrangement. 



to its uncompressed state (from Point C at a = 0.5 MPa 
to D at (7 = 0). The final height of the material after 
the compaction cycle is calculated at Point D, A/io.5MPa> 
and used to obtain the volume fraction at the given stress 
amplitude. The sequence is then repeated for increasing 
values of the mean stress at 1 MPa (D — > E — > F — > G) 
and the correspondingly larger amplitudes of oscillation 
lead to a new height and packing density A/iiMPa for the 
new oscillatory stress. 

The dynamics of compaction is well described by a log- 
arithmic dependence of Ah(t) for each compaction cycle 
valid for long times: 

Ah(t) ln(i). (1) 

At short times, it crosses over to a faster decay. This 
implies that the density behaves as ~ 1/ ln(i) at long 
times in agreement with previous experiments Q- This 
is shown in Fig. \I\&s the solid line from B to C and from 
E to F. The relaxation is so slow that one could argue 
that the final steady-state density has not been achieved 
on the timescale of the experiment. 

The volume fraction after each compaction cycle is 
plotted against the amplitude of oscillatory stress in Fig. 
121 The protocol involves a stepwise ramp of compaction 
cycles from zero stress amplitude to 4 MPa, then back 
to zero and back again to 4 MPa. The figure shows an 
irreversible branch as the mean stress is first increased, 
until a plateau in the volume fraction is reached. As 
the stress is then reduced the system enters into the re- 
versible branch of compaction. By ramping up and down 
in oscillatory stress we find that the system reproduces 
the same curve. These are the reversible jammed states 




Amplitude of oscilation (MPa) 



FIG. 2: Volume fractions obtained by oscillatory compression 
experiments at varying oscillatory stress amplitudes. The fi- 
nal heights at D and G in Fig. are used to calculate the 
volume fractions shown in this figure. 

that we use as the reference states in the stress relaxation 
experiments. 

Stress Relaxation Experiments (strain- controlled). - It 
is important to distinguish between processes related to 
large scale deformations in granular compaction studied 
above and infinitesimal perturbations related to support- 
ing the stress once the volumetric conditions have been 
satisfied. In the latter, an application of an external 
stress will result in dissipation mechanisms quite different 
to the compaction process. Therefore, we next probe the 
mechanisms of energy dissipation of the fully compacti- 
fied system by performing infinitesimal step compression 
experiments and observing the resulting response in the 
stress. 

We perform uniaxial compression tests with the plate- 
cup configuration used above. We contract the system 
by applying a step strain Ae in the range (1 — 3) x 10~ 3 
to the glass bead sample at a given confining pressure. 
Meanwhile, we monitor the temporal evolution of the dif- 
ferential stress, Aa(t). It is defined as the difference be- 
tween the stress at time t measured after the strain is 
applied and the stress before the perturbation is applied. 
It follows that Act(0) is the change in stress between the 
baseline (stress before perturbation) and just after the 
perturbation is applied ||. 

Figurc[3]shows the resulting relaxation at different con- 
fining pressures from 0.1 MPa to 5 MPa. At the slowest 
strain rate of Ae = 6.6 x 10~ 5 s _1 we find a relaxation 
which is logarithmic in time, at long times. When we 
repeat the experiment at the maximum speed allowed by 
the press, Ae = 9.8 x 10~ 3 s _1 , we are able to observe 
the "instantaneous" stress response in addition to the 
subsequent relaxation. In this case we find a two-step re- 
laxation, which is well approximated with the following 
equation, plotted in Fig. [3] for the 4 MPa stress relax- 
ation: 

Acr(t)/Acr(0) = A + Be-^ T1 - Cln(i), (2) 
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FIG. 3: Experiments: stress relaxations at different confining 
pressures. Results are shown for two strain rates. The fastest 
one reveals the two-step relaxation. 



where t\ = 1.4 s is the fast relaxation time and A = 
0.9, B = 0.09 are constants, and C = 2 x 10~ 3 sets 
the rate of the slow relaxation. Since the slow strain rate 
corresponds to a straining time of T = 21.6 s and the fast 
rate to T = 0.18 s for the case of the 4 MPa sample, the 
fast relaxation is only observed in the latter experiment, 
as the relaxation time t± is slower than the application 
of the strain. 

We argue, by analogy with glassy dynamics, that the 
fast relaxation is a single particle relaxation mechanism 
whereas the slow is representative of a collective rear- 
rangement of many grains via sliding and their ensuing 
"aging" properties. We provide supporting evidence for 
this claim by means of molecular dynamic simulations. 

Computer simulations. - In order to decipher the main 
microscopic mechanism, we perform a numerical study 
based on Molecular Dynamics of elasto-frictional spheri- 
cal particles. 

Interparticle forces are computed using the principles 
of contact mechanics and considers: normal Hertz forces, 
F n , tangential Mindlin forces, F t , and dry Coulomb fric- 
tion: F t < [iF n , with [i the friction coefficient. Full de- 
tails are given in |j| . 

We provide two principal mechanisms of dissipation. 
If the grains are touching, they exert contact damping 
forces, which arise from viscoelastic dissipation between 
the grains: f 7 f ss = -7„£ 1/2 £ and ff lss = ~7t£ 1/2 s [HI- 
Reference values for the damping constants 7„ and jt 
can be obtained from [Toj . Furthermore, the grains are 
immersed in a viscous fluid, such as air or water, which 
causes global damping according to the classical Rayleigh 
theory. The drag of a sphere immersed in a viscous fluid 
is F dra9 = 67T7yi?x, where rj is the viscosity of the fluid 
(an analogous expression holds for torque damping) 01 • 

The packings were equilibrated at a given pressure of 
IMPa according to the previously established method Q. 
We probe the macroscopic mechanical properties by ap- 
plying an infinitesimal step strain, and we monitor the 
time dependence of the stress, thus mimicking the cxper- 
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FIG. 4: Simulations: (a) stress relaxation after an instanta- 
neous step strain at IMPa. Red solid triangles correspond 
to a system with global damping n — 1.7 x 10 -3 Pa-s > if, 
blue open circles (also in the inset) to an overdamped sys- 
tem with contact damping 7 n = 2 x 10 -2 kg/s 
2 x 10 -3 kg/s m 1,/2 , and black open squares (also in the inset) 
to a system with critical contact damping. The inset shows 
the behavior of the fraction of the particles sliding in the aging 
regime and at the critical point, (b) Sliding-sliding correlation 
function showing the appearance of aging (dependence on t m ) 
above the critical damping. The black line show the behavior 
below critical damping. Insets show the stretched exponential 
exponent j3 and the characteristic relaxation time tj versus 



imcnts. 

We find a critical damping (both for global and con- 
tact), above which the slow relaxation ensues. The crit- 
ical values are 7^ = 2x 10~ 3 kg/s to 1 / 2 , 7^ — 2 x 10~ 4 
kg/s m 1 / 2 , while the critical global viscosity is rf = 
1.7 x 10~ 4 Pa-s 0. In Fig. 0k we show typical re- 
laxation curves obtained for three systems with cither 
global damping only (77 > ?y c ), or with contact damping 
only (7„ > 7^,7* > 7 t c ), and a third system with criti- 
cal contact damping (7„ =7^,7* = 7 t c ), which dissipates 
energy mainly via Coulomb friction. We see that the 
curves with damping larger than the critical value are 
in good qualitative agreement with those from the ex- 
periment, displaying slow stress relaxation at long times, 
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while the curve with mostly frictional dissipation decays 
very fast. Comparing with real viscoelastic constants for 
common glass materials ^| and the viscosity of air or 
water, we find that real viscoelastic damping is almost 
always above the critical values. Thus, in most experi- 
mental situations, except perhaps in vacuum, the system 
will be ovcrdamped and the slow relaxation will be ob- 
served. 

Considering that we have the full information on the 
motion of the particles, we now investigate the micro- 
scopic origin of the slow relaxation. By monitoring the 
small displacements of the particles we find that the shear 
displacements are key to explaining the relaxation mech- 
anism. The inset of Fig. 0]i shows the fraction of particles 
that are sliding (Ft = fxF n ) in the system at a given time 
(for these particles the shear displacement is significant). 
If the system is undcrdamped (7 < 7 C , for both nor- 
mal and shear), this fraction decays rapidly to zero and 
the stress is quickly relaxed. On the other hand, when 
the system is overdamped (7 > 7 C ) the fraction of slid- 
ing particles increases as a function of time, analogous 
to the strengthening behavior found in previous exper- 
iments 0- This slow dynamical strengthening (which 
may saturate at larger times) is responsible for the slow 
stress relaxation. 

To further investigate the features of the overdamped 
state we study the time correlation function of the sliding 
grains to quantify their motion. At a given time step we 
construct a state-vector s(t) with M components (M is 
the number of contacts in the system); the i-th. compo- 
nent being 1 or -1 according to whether the i-th contact 
is sliding or not. We consider a time correlation func- 
tion C(t,t w ) =< s(t + t w )s(t w ) > with a dependence on 
the waiting time, t w , measured from the time when the 
perturbation is applied. Below the critical damping con- 
dition we find (Fig. that the C(t,t w ) decays rapidly 
with no evidence of t w dependence. With small damping 
forces the grains in the system have no time to develop 
significant shear displacements. The number of sliding 
particles is very small and the system does not display 
any slow relaxation. 

In contrast, above critical damping we find (Fig. 0Jd) a 
stretched exponential decay C(t,t w ) ~ exp[— t/Tf(t w )]P, 
where Tf(t w ) is the characteristic time dependent on the 
waiting time and /3 s» 0.8 is the critical exponent (see 



Insets of Fig. 0Js). Large damping sufficiently slows 
down the system such that the grains have enough time 
to develop substantial tangential displacements, which 
in turn generate shear forces large enough to cause slid- 
ing of the grains. Interestingly, the same effect has also 
been found in the history-dependent behavior of packings 
near the jamming transition, in this case, as a function of 
the compression rate |l3[ . The fact that both damping 
and compression can drive the system into a glassy phase 
through a well-defined transition, which exhibits aging, 
carries important implications. 

Whereas the transition between an undcrdamped and 
overdamped state is not surprising, it is intriguing that 
the latter state show signatures of glassy behavior. Con- 
ventional glassy systems such as polymer melts undergo 
the glass transition by fast cooling of the system, while 
here we show that in granular athermal systems the 
damping plays a similar role to temperature. A granu- 
lar 'glass transition' driven by damping opens interesting 
unifications between the two types of systems, the de- 
tails of which we plan to test experimentally in further 
work. Since damping facilitates dissipation by friction as 
the timescale of grain contacts is prolonged, we predict 
that a packing of frictionlcss droplets in emulsions would 
not experience the observed slow relaxation and aging 
dynamics. In fact our simulations indicate that a system 
with frictionless particles, \i = 0, as well as a system of 
infinitely rough particles, fi — > 00, does not display aging, 
thus confirming that the glassy properties are due to the 
finite friction coefficient of grains. 

In summary, in this work we distinguish between 
stress relaxation processes related to reaching a granular 
jammed equilibrium state and the remaining infinitesimal 
relaxation. Once the system is jammed at a given pres- 
sure, the stress relaxation is characterized by a fast expo- 
nential relaxation followed by a glassy slow logarithmic 
decay. The glassy phase is characterized by the aging of 
the sliding correlation function and the ensuing stretched 
exponential behavior, which arc in turn responsible for 
the slowdown of the dynamics. It's interesting to note 
that the amplification of the sliding could be associated 
with the initiation of an avalanche inside the system. 
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